Exact solutions of the (2+1) Dimensional Dirac equation in a constant magnetic field 

in the presence of a minimal length 
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We study the (2+1) dimensional Dirac equation in an homogeneous magnetic field (relativistic 
Landau problem) within a minimal length, or generalized uncertainty principle -GUP-, scenario. 
We derive exact solutions for a given explicit representation of the GUP and provide expressions 
of the wave functions in the momentum representation. We find that in the minimal length case 
the degeneracy of the states is modified and that there are states that do not exist in the ordinary 
quantum mechanics limit (J3 — > 0). We also discuss the massless case which may find application in 
describing the behavior of charged fermions in new materials like Graphene. 
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I. INTRODUCTION 

In recent years there has been extensive research on 
the minimal length formalism. The concept of a mini- 
mal length has emerged from various studies on quan- 
tum gravity [l|, perturbative string theory and black 
holes Q. See [j] for a recent review. A consequence 
of the presence of a minimal length is that the Heisen- 
berg uncertainty relation becomes modified and this re- 
sults in UV/IR mixing. Consequently it is meaningful 
to study quantum mechanics in the presence of a mini- 
mal length In particular, exact solutions of various 
non relativistic [9l-[l5| and relativistic problems (lll - fl9| 
have been obtained in the presence of a minimal length 
(Axo = hyffi). Approaches have also been discussed that 
try to incorporate a minimal length in the quantum field 
theory formalism and explicit calculations of the Casimir 
effect [13, HH and the Casimir-Poldcr interactions [2l| 
within a generalized uncertainty principle have been de- 
rived. 

Many of the studies available in the literature deal with 
specific calculations, report on the regularizing proper- 
ties of the minimal length, and also have the purpose of 
deriving upper bounds on the minimal length via com- 
parison with experimental measurements, where possi- 
ble. The authors of ref. [23|, [24| for instance solve the 
inverse square potential exactly in arbitrary dimensions 
and show how the minimal length acts as a natural cut- 
off regulator. In ref. [25[ the authors study the scatter- 
ing problem within a GUP scenario and derive the de- 
pendence on hyf]3 of the scattering amplitude and cross 
section. We may note that upper bounds of quite differ- 
ent magnitude have been derived. In ref. [26j the semi- 
classical limit of the GUP scenario has been addressed 
and a quite impressive constraint on the minimal length 
has been derived by computing the perihelion shift in 
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a central force potential. Comparing it with the ob- 
served precession of the perihelion of Mercury results in 
(Ax) m in = tiVfi < 1CU 68 m, some 33 orders magnitude 

below the Planck length (L P = = 1.16 x 10 33 m). 

Other interesting constraints come from including the 
corrections due to the minimal length to the hydrogen 
atom spectrum and computing the Lamb shift. The ac- 
curate measurements available for the Lamb shift allow to 
derive and upper bound on the minimal length of the or- 
der of the electroweak scale: (Ai) mm = h\f]i < 10 -17 m 

[Hill- See also the recent works in and for fur- 
ther discussions about the minimal length phenomenol- 
ogy using a somewhat different GUP representation than 
the one taken up here. 

Here we propose to study a relativistic quantum me- 
chanical problem, namely, the (2 + 1) dimensional Dirac 
equation in the presence of a minimal length. To be more 
specific we shall obtain exact solutions (eigenvalues and 
eigenfunctions) of the Dirac equation in the presence of a 
homogeneous magnetic field (relativistic Landau levels - 
LL-) . This topic has become quite interesting because of 
its application to various branches of physics, particularly 
in condensed matter physics. In passing we may note 
that due to this growth in the interest for 2-dimensional 
electron systems, very recently (non relativistic) Landau 
levels have been for the first time imaged, revealing the 
expected ring-like internal structure of the wave func- 
tions by means of scanning tunneling spectroscopy [3lj |. 
In this context we would like to note that, from the the- 
oretical side, the Pauli equation has also been studied in 
the presence of a minimal length [32]. However, we shall 
implement the minimal length formalism in the first or- 
der Dirac equation rather than after obtaining the second 
order equations for the spinor components. We shall ob- 
tain solutions of the problem after converting the equa- 
tions for the components into Schrodinger like equations 
with some standard solvable potential. Subsequently the 
scalar product for the model (which is quite different from 
the standard one) will be used to determine admissible 
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limits on the angular quantum number m (in the momen- 
tum space) and this in turn will be used to determine the 
spectrum and the corresponding eigenf unctions. A no- 
table feature which emerges from the analysis is that in 
certain cases the admissible values of the angular quan- 
tum number is constrained by a bound which depends on 
the minimal length. Also, there is a class of states which 
cease to exist in the limit f$ — s- 0. Finally it may be noted 
that apart from being interesting in itself, the massless 
Dirac equation in (2 + 1) dimension finds application in 
condensed matter physics. For example, massless Dirac 
equation in (2 + 1) is relevant in describing the motion 
of electrons in graphene [H||. In view of this we shall 
also find the eigenvalue spectrum and the corresponding 
eigenfunctions in the massless case. 

In addition we shall discuss how our results in the 
massless case, relevant to graphene, can be used to ex- 
tract an upper bound on the minimal length by compar- 
ing with experimental measurements of the relativistic 
Landau levels (LL) in graphene as reported in [34| . Our 
upper bound on the minimal length derived from com- 
paring measurements of (electron-electron and electron- 
hole) transitions between the first excited Landau lev- 
els of graphene from [34[ turns out to be (Ax) m i n = 
h\f]5 < 2.3 nm and is of the same order of magni- 
tude of the bound obtained from considerations of the 
corrections due to a minimal length on the Casimir ef- 
fect [2(| • In [35[ the authors use the exact solution of the 
(non relativistic) harmonic oscillator in arbitrary dimen- 
sions within a GUP scenario in order to derive an upper 
bound on the minimal length referring to measurements 
on electrons trapped in strong magnetic fields (Penning- 
trap [36[) whose motion is effectively one dimensional. 
They take advantage of the n 2 dependence of the min- 
imal length correction to the (non relativistic) eigenval- 
ues and derive an upper bound on the minimal length 
(Hy/P < 10~ 16 m) which is however based on the rather 
strong assumption of being able to measure the energy 
eigenvalues up to quite large values of the quantum num- 
ber (n « 10 s ). Their actual bound {h^/j3 < 15xl °~ 9ro ) 
becomes the order of a few nanometers when n ~ 0(1) 
and is of the same order of the bound derived in this work 
(see details in section Ull Ap . 

The organization of the paper is as follows: in scctionim 
we present the problem and obtain the exact solutions; in 
section InTl we analyze the spectrum and provide explicit 
expressions for the momentum space wave functions; fi- 
nally section llVl is devoted to a discussion and conclusion. 



II. (2+1) DIMENSIONAL DIRAC EQUATION 
IN THE PRESENCE OF A MINIMAL LENGTH 
AND WITHIN A CONSTANT MAGNETIC FIELD 

In the minimal length formalism the Hcisenberg alge- 
bra associated with the position coordinates ii and the 
momentum f>i is given by [H @ : 



where j3 > is the minimal length parameter. The corre- 
sponding generalized uncertainty principle (GUP) reads: 

AxiApj > f %[1 + /3(Ap) 2 + /3(p) 2 ] (2) 

yielding a minimal observable length Axq = %\[p. A 
representation of Xi and pi which satisfies Eq.([T]) may be 
taken as 



d 

Xi = ih(l + (3p 2 )—— , pi=Pi 
opt 



from which it also follows that 



AxiAxj > hf3\(piXj —pjXi)\ 
ApiApj > 



(3) 



(4) 
(5) 



It is important to note that the scalar product in this 
case is not not the usual one but is defined as 



(/Is) 



d 2 p 



f*(p)g(p) 



(6) 



Let us now consider the (2+1) dimensional Dirac equa- 
tion in the presence of a homogeneous magnetic field 
B = (0, 0, Bo) with the corresponding Hamiltonian given 
by: 



H = ccr.(P + -A) + <r z Mc 2 



(7) 



where cr = (a x , a y ), and <j z are Pauli matrices and the 
vector potential is chosen in the symmetric gauge: 



Bn 



Bn 



A x — Ay — —x . 

The eigenvalue problem reads: 



(8) 



(9) 



Now using the representation ([3]) the above eigenvalue 
equation can be written as 



/ Mc 2 cP_ 
I cP+ -Mc 2 




[Xi,Pj\ 



ihSij(l 



HP 2 ) 



(1) 



where we have defined 

P±=P X ± iPy = (p x + ± i (p y + (11) 

Written in terms of components, Eq. (fT0| reads 

P_V (2) = P+i> {1) = e+V' (2) , e± = E±M ° 2 

C (12) 

Then decoupling the components we find 

P_P+^ = P+P_V (2) = eV 2) , 

e = (-+(- = -5 (13) 
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Now using the relations ([3]) we find that 



P+ 



j.t> 



P -(i + p P 2 ) [xd p + fd# 

p+(l + /3p 2 ) (\d p -fd 4 



(14) 



where we have defined 

HeBo o 2 2 

A = — — , p x =pcos§, py=psm§, p x +p y =p 

_ (15) 
Following ,5| the generator of rotations in the (x,y) 
plane in the minimal length scenario is defined by: 



xp y - UPx 
1 + Pp 2 



-i hdtf 



(16) 



and satisifies the relations [P±, L z ] = z fhP±. It can then 
be easily verified that the operator J = L z + (h/2)a z com- 
mutes with the Hamiltonian in Eq. (|10[) , so that even in 
the presence of a minimal length we have a conserved to- 
tal angular momentum. Note that in the limit ft — ¥ the 
definition of L z in Eq. (fl~6|) goes smoothly into the ordi- 
nary one. Thus we see that if to is the quantum number 
associated to the operator L z the conserved total angular 
momentum is j = h(m± 1/2). Note that although in this 
instance the angular variable •& is defined in momentum 
space, cf Eq. (|15[) . the quantum number to (associated to 
the eigenfunctions e* of L z ) retains its usual meaning 
of orbital angular momentum quantum number. 

The wave functions may be taken therefore to be eigen- 
states of the (total) angular momentum (note that the 
components have to satisfy the intertwining relations in 
Eq. (|13p ) and we can put them in the form: 



(i) _„(i) 



(pV 



v/ 2 



,(2) 



(p)e' 



i(m+l)iJ 



(17) 



Then from Eq. (|13p we obtain: 



p 2 + 2A (1 + (3p 2 
-A 2 (I + PP 2 ) 2 

p 2 + 2A (1 + I3p 2 
X 2 (1 + /3p 2 ) 



to + 1 — (3 A ( p— to 



dp 2 



ld_ 

p dp 



,21 



x 1 ; 



(p) 



2 ^(P), (18) 



m — /3A p h m + 1 

dp 



dp 2 



1 d 

p dp 



(m + l) 2 



«#(P) = 
(P). (19) 



The above equations are still complicated enough to ad- 
mit direct solutions. However, the solutions may be ob- 
tained readily if we can transform the equations to some 
standard form. To this end we now perform a simultane- 
ous change of wave functions as well as of the variable: 



,(*) 



VP 



■ (0 
tanq, 



1,2 



X ^ [ 2 ' 2 ] 



(20) 



Using the above transformations we obtain from Eq. (|18|) 
and Eq. (QU): 



1 



Ci(C»-l)+6te-l) 



1 r 

2 



dx 2 2 

^te -i)- Cite -i) 



cos 2 (x) 
sin(a;) 
cos 2 (cc) 



where 



1//? 



4/3A 2 ■ 

and the parameters £j and £j are defined as 
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Ci = to + 

C2 = m + 



£1 = m 

6 = to, 



1 

Ha 
1 

rfA 



(21) 
(22) 

(23) 
(24) 



One can identify the above Eq. ([21) as a pair of 
Schrodingcr equations (in units where h 2 /(2M) = 1) with 
the trigonometric Scarf potential of the form: 



V(x) = 



p 2 + v 2 



1 



1 



4 / cos 2 x 



fi 2 — v 2 sin x 
2 cos 2 x 



(25) 



where the parameters \x and v are given in each case 
(i = 1, 2) in terms of the parameters £j and £j via: 



M = & - 



1 



Ci- 



1 



(26) 



We may note that the potential V(x) in Eq. (|25|) has 
certain symmetries that will be of use in writing the so- 
lution of our problem. In particular V(x) is unchanged 
by the replacements /J,,—t — p and/or v — > —v. Upon 
imposing standard boundary conditions on the finite do- 
main x G [— 7r/2, +7t/2] or q € [0, +7r/2] (normalizability 
and vanishing of the wave function at the end-points), 
the eigenfunctions and eigenvalues of Eq. (f/TTj) are read- 
ily obtained from [53, HH : 



ip n (x) = C [z(x)} * [1 - z(x)}%+-z x 

2F1 (—n, /i + i> + l;i/+l;l — z(x)) 



(27) 



1 — sin x 



cos 2 (q) and C is a normaliza- 



K = n 
where z(x) 

tion constant. Note that ipnix) in Eq. (|2"T)) is obtained 
from (37| , via the substitution: p, 1/, z — > 1 — z which 
is easily verified to be a symmetry of the potential V(x) 
in Eq. (25J). 

The vanishing of the wave-function at the end-points 
(i.e. q = and 5 = 7r/2) is ensured by enforcing the 
following constraints: (a) p > —1/2 and v > —1/2; (b) 
fx < 1/2 and v < 1/2; (c) fj, > -1/2 and v < 1/2; (d) 
fi < 1/2 and ^ > —1/2. Solving the parameters p, f 
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m 




1,2 


m > 


(sing) Cl (cosg) Sl 2 F X (-n,n + d +£i>Ci + §;sin 2 g) 


i(2n + Ci+6f 


"I ~TB <™<-l 


(sing) 1 ^ 1 (cosg) Cl 2 Fi (-n,n + l - & + £i,§ - Ci;sin 2 g) 


i ( 2n + l-Ci+c;i) 2 




(sing) 1 ^ 1 (cosg) 1 ^ 1 2 f\ (-n, n + 2 - Ci - £i, | - &! sin 2 g) 


i(2n + 2-Ci-Ci) 2 



TABLE I. ipn.ln and the corresponding energy values for different values of m. In this case a solution to Eq. (|26p is/i = m + l + -^ 
and v = m. The range m > is obtained solving in terms of m the constraints: fi, v > — i (or equivalently Clj£i > 0). The 
range — | — i < m < —1 is obtained solving the constraints: > — |, v < \ (or equivalently (1 — £i) > 0,£i > 0). The range 
m < — 2 — is obtained solving the constraints: < ^ (or equivalently (1 — £i) > 0, (1 — £i) > 0). The fourth possible 
constraint fx < |, f > — | does not admit solutions for any value of m. 



m 


(2) 


fc 2 


m > 


(sing) C2 (cosg)« 2 2 Fj (-n, n + £ 2 + 6, C2 + §;sin 2 g) 


1 (2n + C2 + 6) 2 


-| - -ts < m < -1 


(sing) W2 (cosg) ?2 2 Fi (-n, n + 1 - Ca + €2, | - C2;sin 2 g) 


1 ( 2n + l-C2+6) 2 


m <\-TB 


(sing) 1 ^ 2 (cosg) 1 ^ 2 2 f\ (-n,n + 2 - (2 - £2, § - C2; sin 2 g) 


i (2n + 2 - C2 - 6) 2 



TABLE II. ip„, m and the corresponding energy values for different values of m. In this case a solution to Eq. (|26p is /i = m+ 4r; 
and 1/ = m + 1. The three ranges of the m values are found solving the same constraints described in the caption of TableHJ 
Note that in this case the fourth constraint p, < |, v > — | has solutions for m in the range — § < m < 5 — ^ wmcn is 
meaningful only for ^ < 2. However the minimal length is physically expected to be a small quantity and we have indeed 
-ij » 1. See the discussion in the text and Eq. (|35[) . So this possibility will be ignored throughout. 



in Eqs. (I26|) in terms of the angular momentum quantum 
number m provides with the three ranges (of m) in Tables 
U and HIl Note that one of the constraints does not have 
solution for any value of m. The wave-functions in Tables 
HI and [Til are then obtained from Eq. ((27)) . In the second 
and third row of both tables repeated use is made of the 
fact that the potentials in Eq. (|21[) are invariant under 
the reparametrization: 



C -> 1 - Ci , and/or & -> 1 - & 



(28) 



We conclude this section with a final important re- 
mark. While we have applied standard boundary condi- 
tions in the finite x (or q) domain (normalizability and 
vanishing of the wave function at the end-points) it is 
interesting to note that these can be transported back to 
the physical (radial) p-spacc of our original Dirac prob- 
lem and can be given a physically sound interpretation. 

The normalization integral of the Dirac spinor reads: 



d?p 

l + (3p 2 



+ U {2) ) v 



',(2) 



(29) 



and normalizability of the spinor solution is guaranteed, 
if both radial components components satisfy: 



pdp 



'(p)| S 



< 00 



(30) 



Because of the deformation of the measure, introduced 
by the minimal length, the asymptotic behavior of the 
u^ l '(p) functions that ensures such condition is: 



u«(p) 



1 

00 p x 



X>0 



(31) 



In the reduced problem, where ip^ (p) = p 



less x is large enough (x > 1 /2) the wave function ip^ (p) 
will not vanish as p — > 00. Thus in this sense we con- 
clude that normalizability alone of the (p) wave func- 
tions does not warrant that the reduced wave-functions 
<p( l >(p) vanish at p — > 00 (or in g-space at q — 7r/2), 
while (p^(q)\ _ w = is the standard boundary con- 

I q— 2 

dition which we have implemented in building up the 
results of Tables HI and HIl Note that for a vanishing mini- 
mal length (/? — > 0) the measure reduces to the standard 
one and the normalizability of the (p) requires instead 
X > 1 which would ensure that the reduced wave function 
vanishes as p — > 00 (or q = ir/2). 

On the other hand in our relativistic Dirac problem 
the energy integral computed from the quantum Dirac 
hamiltonian in Eq. (fTCTj) is: 



d 2 p 



Mc 2 (V^V 1 ^ 



l + /3p 2 

(V 1 ))* cP.^-Mc 2 (V 2 ))V 2 ) + 



(> 2 >) cP+V> (1) 



(32) 



Now require in addition the finiteness of the energy in- 
tegral. This means that the second and fourth integrals 
in Eq. (|3"2"1) must be finite. The operators P± as given in 
Eq. (|14[) are linear in the radial momentum p (as it is ex- 
pected from a Dirac Hamiltonian). Then assuming that 



as p — > 00 the components behave as u 



p x (with 



X > to ensure normalizability) the asymptotic beahv- 
ior of the integrands in the second and fourth integral in 
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Eq. (02) is: 



pp z p x p x 



which turn out to be integrable only if 2\ > 1 or \ > 1/2 
which, for example in the case of the first row of Table U 
(X = £i + 1/2), translates into £i > or /i > —1/2 
ensuring that tp^l _ 2L =0. 

Similar considerations can be performed as regards the 
behavior of the wave- functions at p — 0, and again the 
vanishing of the <pW (q) | is ensured by the finiteness 
of the energy integral. Note that in the limit p — > the 
operators P± in Eq. [TJ] will be dominated by the deriva- 
tive terms (d p ) which give an extra inverse power of p in 
the third and fourth integrals of Eq. 1321 In particular we 
have verified that all the conditions discussed to deduce 
Tables U and [TT] can be deduced in the radial p-space by 
imposing the finiteness of the energy integral. 

Let us conclude these considerations with a final obser- 
vation. Physically we would expect that the wave func- 
tions u^'(p) and thus the ip^(p) have always a regular 
behavior at p = 0. However this is not excluded by the 
boundary condition that we have imposed in the reduced 
problem. From Eq. (fMf we see that the condition to 
require that the function would not be divergent at 
p = is d — 1/2 > while the condition that we have im- 
posed is the less restrictive one > 0. We notice however 
in Tables Hill IIVI and IVl that our wave functions never di- 
verge for p — > 0. This can be understood by the fact that 
the Q of our problem are not continuous parameters but 
are instead discrete because they depend on the orbital 
angular momentum quantum number m. Indeed in the 
case of Eq. [34] ( = m + 1/2 and the condition £ > 
reduces to m > —1/2 which is effectively equivalent to 
m > (since m is integer) or £ — 1/2 > 0. 



III. SPECTRUM AND WAVE FUNCTIONS IN 
MOMENTUM SPACE 



Since we have reduced ourselves to the exact study of a 
Schrodinger equation in the x (or q) space, the problem 
does not need any further inspection. However, for a 
better understanding as well as for completeness the full 
spinorial solutions will be given in the p-space. Starting 
from tables I and II and the wave functions in the 
q space, we simply obtain the form of the corresponding 
"radial" wave functions in the p-space through Eq. (fl7|) 
and Eq. (f2"0)) . They arc found to be of the form (C, is a 
normalization constant) 



,(P) = c t 



(i+/V) 



2 -Pi 



-n, Ci + & + n, Ci + 



i /V 

2' 1 + (3p 2 



(34) 



where Q and are defined by Eq. (1231) and Eq. 
and 2-Fi is the hypergeometric series [39|]. We shall now 
classify the eigenfunctions and the corresponding energy 
values according to the angular quantum number m. The 
results are summarized in tables IIII1 IIVI and [V] We give 
here explicitly only the positive branch of the spectrum. 
For the negative eigenvalues similar formulas are readily 
obtained (see Eq. (|2"2"])). 

Then, we build up the full spinor solutions by putting 
together those states of V'*' 1 ' an d ip^ which have the 
same energy eigenvalues. In doing this one has to use the 
last columns of Tables HI and [TT1 It turns out that the first 
row of Table U provides the same energy spectrum as the 
corresponding row of Table ITO together, they produce Ta- 
ble IIIII Analogously for the second and third rows which 
respectively imply Tables IIVI and [V] One can verify that 
the upper and lower spinorial components are actually 
interlaced by (|12j) . The ranges for the quantum number 
m in Tables IIIII IIVI and [V] descend from the first column 
of Tables |T] and [TTJ When coupling the first row of Table 
U and the first of Table [H] one takes the intersection set 
of the two m-domains. The same is done for the other 
rows. This produces the three different possible sets of 
values of m presented in Tables IIIII IIVI and [V] 

Note also that in Eq. (O and tables UTTT HVl and M 
the normalization constants Ci, with (i — 1,2), are not 
independent. This can be seen from the fact that the 
two spinor components ip^ must satisfy the intertwin- 
ing relations (|12p or from the fact that requiring a given 
normalization for the full spinor solution gives a relation 
between C\ and C%. 

Let us now examine and discuss the results presented 
in these Tables. 

Table Hill gives the energies and the eigenfunctions with 
positive values of the angular momentum quantum num- 
ber (m > 0). We find that all the energy levels except the 
lowest state (which is a singlet) have a finite degeneracy. 
For example, for n + m = N the levels are (N + l)-fold 
degenerate. Also, all states are doublets i.e. have a spin 
up as well as a corresponding spin down component. Ta- 
ble |IV] shows the results for values of the angular momen- 
tum quantum number in the range 

Note that: 



.1 __L.< m <_i. 



2c 



2hMc 



1 

I 



= 2 



lr. 



Ax 



> 1, lc 



Mlol 



(35) 



where lol = \f - denotes the electron cyclotron frequency 
and l c is just the characteristic length of the associated 
oscillator, which has to be considerably larger than the 
minimal observable length if this very problem has to 
be studied. Therefore m can not assume an arbitrar- 
ily low negative value but is constrained by the lower 
limit (— i — jp) (which is a very large negative num- 
ber by virtue of (|35|) and in the limit j3 — > it becomes 
infinitely negative). These class of energy levels have a 



6 



TABLE III. Energy levels and the corresponding wave functions for m > 0. A given energy level with n + m — N has a finite 
degeneracy D = N + 1. 



i[l +J gSga( n + m )] n = 0,1,2,... ip n , m = ( ^' 



P m o-Fl I -w.n+2(m+l)+-ra,m+l, |9p2 ., I , . „( m + 1 ') p. ( _ n n ^2(?n+l)+ -A. ,m+2, ^ p , 2 ,. | 



ft.m = Cl — l -e = C2 



TABLE IV. Energy levels and the corresponding wave functions for — ~ — jjj < m < — 1. The degeneracy D of these levels is 
finite and explicitly given by D = + -^]. In the limit of a vanishing minimal length (/3 — > 0) the degeneracy of these energy 
levels becomes infinite (D — ^ 00). 



£0 = Mc 2 



00, 







-En 



M 2 c 4 + 2heB cn (1 + /3 2 
n = 1,2, 



^n,m - \ ^( 2) 



Fl (_„,„+ 1 ~Z7ZSZ] 

2 1^ ■ ''l + /3p* J j( m+ l)fl 



2 F 1 (- w ,Ti+2+^,|m|+l,- 



iij- 



1 



(2) ^ 
V'n.m = C 2 - 



TABLE V. Energy levels and the corresponding wave functions for m < — | — In this case, similarly to what happens 
for the levels in Table [TTT1 the degeneracy of the energy levels with n + |m| = N and AT > [~ + 3^5] + 1 is finite and given by 



D = N-[i + ^) 



E n , m = JaPc* + 2heB c(n + \m\) [/3^ (n + \m\) - l] n = 0, 1, 2, . . . = ( ^p™ J 



P 1 "' 1 9-F1 f -n,n+2|m|- A,l+|m|, ,3p T | „i|m|-l) F / +2 i i_ l i i gp 

fp'n , m — ^ 1 ] 1 ~ 1 1 G (^71 , m — "—'2 " j ^ j 1 c 



(l+/3p2)' m| W (l +| 9„! 



finite degeneracy £> = [i + A-] for finite values of /?. 
£> becomes infinitely large when (3 — > and this family 
of states reduces to the ordinary quantum states of the 
relativistic Landau problem with negative values of m. 
Interestingly the ground state is a spin singlet while the 
excited states are spin doublets. Table [V] gives the en- 
ergy eigenvalues and cigenfunctions with m in the range 
m < — \ — jp. We note that this range becomes mean- 
ingless when (3 — > and the corresponding states loose 
therefore any physical meaning in this limit. However for 
finite values of /3 (a non zero minimal length) such states 
are physical states and must be included in the physical 
spectrum. They are all doublet states and the energy 
levels for which n + \m\ = N with TV > [± + -X] + 1 the 

degeneracy is given by: D = N — [i + j^]. 

From the above tables it can also be seen that for 
m = — \ — there isn't any acceptable spinorial so- 
lution. This directly descends from our g-space analysis 
where one observes that such a value of m would make 



it necessary to appeal to the second line in Table Q] (for 
the upper component) and to the third line in Table ITT1 
(for the lower one). These solutions cannot be coupled 
though, as it is straightforward to verify that they don't 
share the same energy, or in other words that the corre- 
sponding p-space components and tp^ thus obtained 
do not verify p^|) . 

We observe here that the no-GUP context is correctly 
reproduced by letting j3 — > 0, because in this case Ta- 
ble |V] along with its angular momentum domain of va- 
lidity becomes meaningless, and the degeneracy D for 
the negativc-m solutions of Table IIVI approaches infinity. 
Note also that in the mentioned limit the following rela- 
tion holds between the hypcrgcomctric and the confluent 
hypcrgcometric series 

(36) 

as can be straightforwardly checked from their standard 
definitions 



2Fl (a,b,c;z) = ^-^—- lFl (a,c; X ) = ^—- \z\ < 1 (37) 

fc=0 v ' fc=o v 1 
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so that 

and the expected eigenfunctions of the ordinary quantum 
mechanical treatment are obtained. 

A. Massless Dirac equation in (2 + 1) dimensions 

The equation governing the motion of electrons in 
graphene is similar to the (2 + 1) dimensional massless 
Dirac equation except that the electrons move with Fermi 
velocity vf ~ 3^0 c rather than with the velocity of light 
c. The Hamiltonian for the electrons in graphene in the 
presence of a magnetic field is 

H D =v F (T-{p+-A) (39) 
c 

In order to obtain the spectrum and the wave functions 
one needs to make minor changes in the results obtained 
earlier. The energy levels are presented in Table IVII (we 
only write down the spectra since all the rest stands un- 
changed). 

From Table |VI] we find that for to > 0,n + m = N 
the levels are (N + 1) fold degenerate. For —\ — < 
m < — 1 there is a zero energy ground state which is a 
spin singlet while the excited states are spin doublets. All 
these states have a (large) finite degeneracy with respect 
to m which is as before given by D = [5 + j-y]. Finally, 
for to < — I — -Aj the degeneracy of the levels with n + 
\m\ = N, N > ^ + -X is similarly calculated to be given 
by: D = N — [| + -^). As in the massive case these 
solutions become meaningless in the limit of a vanishing 
minimal length (f3 — > 0). Finally to see how the energy 
levels deviate from the usual relativistic Landau levels, 
let us examine the spectrum in the small /3 limit. For 
this let us consider the second energy level in Table VI 
and expanding with respect to j3 we find 

1 2heB n heB Q n I heB n 

E ^ = VF \j—^ n+l3vF ^ry^- + °(P ) ( 4 °) 

where the first term gives the usual Landau levels. Thus 
the spectrum contains an additional term involving n 3 / 2 
and the dispersion relation is indeed modified. These 
modifications of the dispersion relation might have im- 
portant implications, for instance in the calculation of 
quantities like the density of states, which will be ad- 
dressed elsewhere. Here we would like instead to com- 
ment on how our exact solution of the GUP Dirac equa- 
tion in 2+1 dimensions in a constant magnetic field might 
be already of use in deriving an upper bound on the min- 
imal length. 

In order to do so we can compare our derived formulae 
for the energy spectrum with the experimental measure- 
ment of the transitions between graphene Landau Levels 
(LL) H, 



In ref. [H| infrared spectroscopy has been used to re- 
solve the transitions between graphene LL in the pres- 
ence of magnetic fields of intensities up to Bq = 18 Tesla. 
The authors report a linear behavior of these transmis- 
sion resonances with ^/Bq from which a best fit value 
of the fermi velocity v F = (1.12 ± 0.02) x 10 6 m/s is 
deduced. From this experimental value of the Fermi ve- 
locity one can deduce for instance for the first excited 
level of the graphene Landau spectrum in the absence of 
a minimal length (c.f. first term of Eq. (j40|) with n=0), 

E = VfJ^^- for B = 18 T the (experimental) value: 
E = (172 ± 3) meV (41) 

and note that (SE)/E = (5vf)/vf- We wish to use this 
result to provide an upper bound to the observable min- 
imal length Axo — hy/]3 of our GUP model. We may use 
the results of Table VI. Picking the spectrum on the first 
line and setting to = 0, n = 1 one has for the energy of 
the first graphene LL: 

The impossibility to experimentally distinguish the devi- 
ation brought about by the existence of a minimal length 
means that the two values, predicted (e[^q) and exper- 
imental {Ef^ ~ E) must be close enough, i.e. they 
must be, with respect to each other, within the experi- 
mental error (c.f. Eq. (|41|Q hence we can surely assume 
that: 

AE = E { $ ~ E ( l= tt) < 6 meV (43) 

from which: 

= E (y/l + S - l) < 6 meV (44) 

where we have defined S = jS^ffa = {hy/fi) 2 |fa. Hence 
since 5 is expected to be a very small quantity we obtain 
the constraint: 

& < ^ « 0.07 (45) 

which in turn, resorting to gaussian units, leads us (with 
B = 18 T = 18 x 10 4 Gauss) to: 

Ax = Ht//3 < 2.3 nm (46) 

thus providing in principle an upper bound on the mini- 
mal length (or equivalently on the parameter j3) appear- 
ing in the framework of a generalized uncertainty princi- 
ple. 
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TABLE VI. Energy levels for massless electrons. The degeneracy of the the energy levels in the massless case is similarly 
discussed as in Tables IIII1 IIVI and [V] . 





m > 


n 


= 0,1,2,... 






(n + m) [i +J gaga(„ + m )] 


-^-■^<m<-l 


n 


= 


E = 







n 


= 1,2,... 


E n = 




;i + /?^n) 


m < ~5 - X/3 


n 


= 0,1,2,... 


E n , m 




(n+lmD^^Cn+H)-!] 



As a final remark we wish to point out that in ref. [34| 
the authors find some discrepancies on the value of the 
Fermi velocity deduced from different LL transitions and 
warn about possible difficulties of the simple interpreta- 
tion of IR data in terms of a simple LL energy subtraction 
based on standard one particle quantum mechanical re- 
sults and conclude that many particles effects may be ex- 
pected to contribute to the LL transition energies. These 
many particle effects may therefore also affect the upper 
bound derived here (c.f. Eq. (|46[) ) on the minimal length. 

Admittedly the upper bound in Eq. (|4"o]l is not a very 
strong bound. It is however comparable with those de- 
rived in [20| where the upper bound obtained from the 
Casimir effect for the minimal distance of the plates of 0.5 
Jim, depending on the particular GUP model, is in the 
range ~ 29 — 58 nm. Similar order of magnitude upper 
bounds on the scale of p hase-space non commutativity 
have been derived in pOj . 

We could perhaps note that a possibility to make our 
bound more stringent would be to follow the approach of 
ref. [35| and assume that in future experiments it will be 
possible to measure LL transitions for very large values 
of the quantum number n. Our argument that led to 
Eq. (|4"o]l could be reproduced for the n-th LL level and 
would provide the bound: 

y n (E n ) cxp eB Q y/n 

(assuming, somewhat optimistically, the same value for 
the relative error of the measure of E n as for the first 
excited level, c.f. (|45|) ). As discussed in Ref. [4l[ a rather 
constraining upper bound on the minimal length comes 
also from the hydrogen IS — 2S transition: ( Ax) m j n = 
hy/]3 < 1CU 2 fm =1CT 17 m. The authors of [Ilf argue 
that this bound could be avoided by assuming that the 
parameter f3 is not a universal constant and could vary 
from one system to another depending, for example, on 
the energy content of the system (the mass of the par- 
ticle, for instance) or the strength of some interaction. 
Indeed, in [411 ] . by making this hypothesis, the authors, 
through a comparison with the experimental results for 
ultracold neutron energy levels in a g ravitational quan- 
tum well (GRANIT experiment) |42|, derive a relaxed 
upper bound to the minimal length which turns out to be 
of the order of a few nanometers [(Ai) mm < 2.41 nm], 
which is quite close to the one derived here (c.f. (1431) ). 



Clearly we could as well advocate the non universality 
of j3 in order to evade the stronger constraints as those 
discussed in and also in [13, HU . 



IV. DISCUSSION AND CONCLUSION 

We have obtained exact solutions of the (2 + 1) dimen- 
sional Dirac equation in an external homogeneous mag- 
netic field in the presence of a minimal length. We work 
within a momentum space representation of the Heisen- 
berg algebra and through an appropriate transformation 
of both the wave function and the variable the second or- 
der equations for the Dirac components are reduced into 
a finite domain Schrddingcr like exactly solvable problem 
(trigonometric Scarf potential). Interestingly it is shown 
that the ordinary boundary conditions in the finite do- 
main (vanishing of the wave function at the end-points) 
can be transported back to the radial p-spacc and inter- 
preted in terms of the finiteness of the energy integral. 

The solutions show that a non-zero minimal length 
changes the spectrum to a large extent as compared to 
the standard relativistic Landau problem. A notable fea- 
ture of this problem is that when the angular momentum 
quantum number m is negative it is constrained and dif- 
ferent ranges of its value point to different class of physi- 
cal states. However, the constraint on the quantum num- 
ber m disappears as the minimal length vanishes (/3 — > 0). 
This can be seen in Tables IIVI and [V] Another feature 
worth noting is the degeneracy pattern of the energy lev- 
els. In the usual relativistic Landau problem, the Landau 
levels are infinitely degenerate for m < 0. In contrast, in 
the present case some energy levels are finitely degener- 
ate (as in Table IIII[) while others have very large finite 
degeneracy (Table IIV[) . It may also be noted that an in- 
teresting feature of the minimal length scenario turns out 
to be the appearance of the solutions reported in Table 
IVl Indeed these solutions exist only for (3^0. In the 
limit (3 — >• the related range of m becomes meaningless 
and also the corresponding eigenfunctions are no longer 
physically acceptable. In this limit, the correct non min- 
imal length situation can be recovered from the results 
of Tables |m] and IVl 

Wc have briefly discussed how our exact solution of the 
problem in the massless case might be used to provide 
an upper bound on the minimal length via a compari- 
son with exisiting experimental measurements of transi- 
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tions between graphenc LL. Finally we wish to point out 
that, always in the massless case, it would be of interest 
to compute other physical quantities e.g, Hall conductiv- 

I 



ity where a comparison with experimental results may 
provide perhaps more stringent bounds on the minimal 
length hyf]3 than those discussed here. 
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